home *** CD-ROM | disk | FTP | other *** search
/ Mac-Source 1994 July / Mac-Source_July_1994.iso / C and C++ / Libraries / stringsearch / bmsource / no.rev.d12.c < prev    next >
Encoding:
C/C++ Source or Header  |  1991-05-06  |  2.7 KB  |  138 lines  |  [TEXT/MPS ]
  1. /*
  2.     search routine generated by gen.
  3.     skip=no, match=rev, shift=d12
  4. */
  5. /*
  6.  * The authors of this software are Andrew Hume and Daniel Sunday.
  7.  * 
  8.  * Copyright (c) 1991 by AT&T and Daniel Sunday.
  9.  * 
  10.  * Permission to use, copy, modify, and distribute this software for any
  11.  * purpose without fee is hereby granted, provided that this entire notice
  12.  * is included in all copies of any software which is or includes a copy
  13.  * or modification of this software and in all copies of the supporting
  14.  * documentation for such software.
  15.  * 
  16.  * THIS SOFTWARE IS BEING PROVIDED "AS IS", WITHOUT ANY EXPRESS OR IMPLIED
  17.  * WARRANTY.  IN PARTICULAR, NEITHER THE AUTHORS NOR AT&T MAKE ANY
  18.  * REPRESENTATION OR WARRANTY OF ANY KIND CONCERNING THE MERCHANTABILITY
  19.  * OF THIS SOFTWARE OR ITS FITNESS FOR ANY PARTICULAR PURPOSE.
  20.  */
  21.  
  22. #ifndef    CHARTYPE
  23. #define    CHARTYPE    unsigned char
  24. #endif
  25. #define    MAXPAT    256
  26.  
  27. #include    "stats.h"
  28.  
  29. #ifndef    TABTYPE
  30. #define    TABTYPE    long
  31. #endif
  32. typedef TABTYPE Tab;
  33.  
  34. static struct
  35. {
  36.     int patlen;
  37.     CHARTYPE pat[MAXPAT];
  38.     Tab delta1[256];
  39.     Tab delta2[257];
  40. } pat;
  41.  
  42. prep(base, m)
  43.     CHARTYPE *base;
  44.     register m;
  45. {
  46.     CHARTYPE *skipc;
  47.     register Tab *d2;
  48.     register q1, tp, t, qp, jp, kp;
  49.     Tab f[256], f1[256];
  50.  
  51.     pat.patlen = m;
  52.     if(m > MAXPAT)
  53.         abort();
  54.     memcpy(pat.pat, base, m);
  55.     skipc = 0;
  56.     stats.len = m;
  57.     d2 = pat.delta1;
  58.     for(j = 0; j < 256; j++)
  59.         d2[j] = m;
  60.     for(j = 0; j < m; j++)
  61.         d2[base[j]] = m-1-j;
  62.     d2 = pat.delta2;
  63.     for(j = 1; j < m; j++)
  64.         d2[j] = 2*m-j;
  65.     for(j = m, t = m+1; j > 0; j--, t--){
  66.         f[j] = t;
  67.         while((t <= m) && (base[t-1] != base[j-1])){
  68.             if((m-j) < d2[t])
  69.                 d2[t] = m-j;
  70.             t = f[t];
  71.         }
  72.     }
  73.     q1 = t;
  74.     t = m+1-q1;
  75.     qp = 1;
  76.     for(jp = 1, kp = 0; kp < t; jp++, kp++){
  77.         f1[jp] = kp;
  78.         while((kp >= 1) && (base[jp-1] != base[kp-1]))
  79.             kp = f1[kp];
  80.     }
  81.     while(q1 < m){
  82.         for(j = qp; j <= q1; j++)
  83.             if(m+q1-j < d2[j])
  84.                 d2[j] = m+q1-j;
  85.         qp = q1+1;
  86.         q1 += t-f1[t];
  87.         t = f1[t];
  88.     }
  89. /*for(j=1; j<=m; j++)printf("[%d]=%d ", j, d2[j]); printf("\n");/**/
  90.     d2[0] = m+1;        /* the case where the match succeeded */
  91. }
  92.  
  93. exec(base, n)
  94.     CHARTYPE *base;
  95. {
  96.     int nmatch = 0;
  97.     register CHARTYPE *e, *s;
  98.     register int s_offset;
  99.     register CHARTYPE *p, *q;
  100.     register n1 = pat.patlen-1;
  101.     register CHARTYPE *ep;
  102.     register Tab *d1 = pat.delta1;
  103.     register Tab *d2 = pat.delta2+1;
  104.     register k1, k2;
  105.  
  106.     s = base+pat.patlen-1;
  107.     e = base+n;
  108.     s_offset = 1-pat.patlen;
  109.     ep = pat.pat;
  110.     while(s < e){
  111. #ifdef    STATS
  112.         stats.slow++;
  113. #endif
  114. #define    RH    (s+n1+s_offset)
  115.         for(p = pat.pat+n1, q = RH; p >= ep; ){
  116. #ifdef    STATS
  117.             stats.cmp++;
  118. #endif
  119.             if(*q-- != *p--){
  120.                 q++, p++;
  121.                 goto mismatch;
  122.             }
  123.         }
  124.         nmatch++;
  125.     mismatch:
  126.         k2 = d2[p-pat.pat];
  127.         k1 = d1[*q];
  128.         if(k2 < k1)
  129.             k2 = k1;
  130.         k2 = q+k2-RH;
  131. #ifdef    STATS
  132.         stats.step[k2]++; stats.jump++;
  133. #endif
  134.         s += k2;
  135.     }
  136.     return(nmatch);
  137. }
  138.